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Abstract 

The main purpose of this paper is to discuss the link between forward-backward 
multiplicity correlations properties and the shape of the corresponding final charged 
particle multiplicity distribution in various classes of events in different collisions. 
It is shown that the same mechanism which explains the shoulder effect and the H n 
vs. n oscillations in charged particle multiplicity distributions, i.e., the weighted su- 
perposition of different classes of events with negative binomial properties, repro- 
duces within experimental errors also the forward-backward multiplicity correla- 
tion strength in e + e~ annihilation at LEP energy and allows interesting predictions 
for pp collisions in the TeV energy region, to be tested at LHC, for instance with 
the ALICE detector. We limit ourselves at present to study substructures properties 
in hadron-hadron collisions and e + e~ annihilation; they are examined as ancillary 
examples in the conviction that their understanding might be relevant also in other 
more complex cases. 
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1 Forward-backward correlations in e + e and hh collisions. 

Forward-backward (FB) multiplicity correlations have been studied in hadron-hadron collisions 
and e + e~ annihilation [Q- 0' S- B @> © H> © 01 - ^ nas ^ een f° un d mat these correlations are much 
stronger in hadron-hadron collisions than in e + e~ annihilation. 

Let us start with UA5 Collaboration results at CERN pp Collider at 546 GeV cm. energy 
[jj], 0] . It was found, by studying in each event the number of charged particles falling in the 
forward hemisphere, n F , and in the backward hemisphere, n B , that the relation between the 
average number of charged particles in the backward hemisphere, n B , and n F is very well 
approximated by a linear one, 

n B (n F ) = a, + bn F , (1) 

where a is the intercept on the vertical axis and b the slope of the linear fit, i.e., the correlation 
strength: 

u_ {{n F - n F ){n B - n B )) 



[((n F -n F y)((n B -n B y)] l/2 ' 

Here the forward hemisphere corresponds to the region of the outgoing proton and the backward 
hemisphere is the symmetric region in the opposite direction. In order to avoid correlations due 
to kinematical constraints like phase- space limits and energy momentum conservation occur- 
ring at the border of the rapidity range available in the collision, and short range correlations 
produced from particle sources in more central rapidity intervals, the study has been performed 
in the pseudo-rapidity interval 1 < \r]\ < 4. It has been found that b parameter is equal to 
0.43 ± 0.01, a much larger value than that found at lower energy (ISR and EHS) [0, 0, |5p, e.g., 
b = 0.156 ± 0.013 at 63 GeV cm. energy. Assuming uncorrected random emission of charged 
particles in the selected intervals of the pseudo-rapidity axis, the shape of the n F multiplicity 
distribution at fixed full multiplicity n = n B +n F is binomial with probability p = 1/2 for a par- 
ticle of the full sample, n, to fall in the backward or forward hemisphere, and thus has variance 
d\ F (n) = p(l — p)n = n/A. Since the experimental value of the dispersion of the full distribu- 
tion, D n , in the above mentioned pseudo-rapidity interval is 9.2 ± 0.1, and n = 15.8 ± 0.2, one 
has 



£S-4<d*») Dl-n 



„ _ npy " = —H — - = o 69 (3) 

Dl + A(dl F {n)) Dl + n ■ 

a much higher value than the experimental one (b = 0.43 ± 0.01). This fact led the Collab- 
oration to assume that not particles but particle clusters of approximately the same size, M, 
are binomially distributed in the two hemispheres and that the decay products of each cluster 
remain within the same hemisphere. Accordingly 

4<4») = Mn. (4) 

But a reasonable agreement with experimental data is obtained by calculating M eS for clusters 
of a mixture of sizes. One finds that 

, D 2 Jn-M eS , n2 



Dl/n + M eff 



with M eff = M cluster + Di ster /M cluster . (5) 
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Here M c iuster and -Duster 316 respectively the average charged multiplicity within a cluster and 
its dispersion. An even better agreement is obtained by allowing a certain amount of particle 
leakage from one hemisphere to the other. In conclusion the correlation strength increases with 
energy and it is the result of binomially distributed clusters, of approximately 2.2 particles per 
cluster, in the two hemispheres. 

Forward-backward correlations have been studied also by OPAL Collaboration in two- 
jet and three-jet events in e + e annihilation at LEP. The forward hemisphere is here chosen 
randomly between the two defined by the plane perpendicular to the thrust axis; although this 
definition is a bit misleading, a smaller effect than that seen by UA5 collaboration is observed. 
FB multiplicity correlations are absent in the two separate samples of events but it turns out that 
a correlation strength b = 0.103 ± 0.007 is found when they are superimposed, a result com- 
parable with that found by the DELPHI Collaboration (b = 0.118 ± 0.009) [JTp. By comparing 
these results with that found at TASSO [@], e.g., 0.084 ± 0.016 at 22 GeV, we can conclude that 
in e + e~ annihilation the energy dependence of the correlation strength is quite weak. 

Our study is motivated by the just mentioned experimental facts and by the finding that the 
superposition of weighted negative binomial (Pascal) multiplicity distributions (NB (Pascal) 
MD's), each describing a different class of events, (soft and semi-hard events in pp collisions, 
2- and 3-jet events in e + e~ annihilation) explains quite well the characteristic features of global 
event properties of collisions in the GeV region, like the shoulder effect in the total n-charged 
multiplicity distributions, P n , and the H n vs. n oscillations (H n is here the ratios of n-particle 
factorial moments, F n , to the n-particle factorial cumulant moments, K n ). The main purpose 
of this search is to show that FB multiplicity correlations in pp collisions and in e + e~ annihi- 
lation as well as their different behaviour can be understood in terms of the same cause which 
explained P n vs. n and H n vs. n general properties, i.e., the superposition of different sub- 
structures. In addition, the assumption that each substructure (class of events) is described by a 
NB(Pascal)MD allows sound quantitative predictions on parameters which are not known from 
experimental data. Coming to the correlation strength obtained in Eq. (]3|), for instance, if par- 
ticles are independently produced and binomially distributed in the two hemispheres, and the 
overall MD is a negative binomial with characteristic parameters n, the average charged multi- 
plicity, and k (it is linked to the dispersion D n by the relation (D 2 — n)/n 2 = 1/k), one gets 
this simple Equation 

^ n 2 /k + n — n n 
n 2 / k + n + n n + 2k 

A formula which can be applied also to individual substructures satisfying the requested condi- 
tions. 

The use of the NB(Pascal)MD in high energy phenomenology is not indeed an arbitrary 
artifact but a consequence of our ignorance on a problem (lack of experimental data and of 
explicit QCD calculations) and is quite often supported by fits with excellent chi squares for 
charged particles multiplicity distributions of the substructures appearing in the total charged 
particle multiplicity distributions through the above mentioned anomalies. 

The occurrence of the NB(Pascal)MD distribution in high energy phenomenology has been 
interpreted since long time as the by now quite well accepted idea that the production process 
is a two-step process: to an initial phase in which a certain number of sources, which have 
been called clans, are independently emitted, it follows their decay into final particles. All 
correlations among produced particles originated by the same source are exhausted within the 
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same clan. The average number of clans N is related to the standard NB parameters, n and k, 
by the following relation 

N = Hn(l +fi/k). (7) 

Accordingly 

D 2 n = nexp(N/k). (8) 

Each clan is either forward or backward. Produced particles by each clan stay all in the same 
hemisphere where the clan is or some of them may leak to the opposite hemisphere: the no- 
leakage and leakage cases will be discussed in Sec. ^] under the assumption that substructures 
in hadron-hadron collisions are described by NB(Pascal)MD and therefore the concept of clan 
rather than the concept of cluster should be used. 



2 Superposition of different classes of events and its influence on the correlation strength 

In this Section the effect of the superposition of different classes of events and its consequences 
on the strength b of FB multiplicity correlations will be discussed independently from the as- 
sumption that charged particle MD's of the two classes of events are of NB type. Results are 
therefore of general validity and can be applied to any pair of MD's describing experimental 
data in the two substructures of the total MD. Accordingly in the following the number 1 and 
2 will indicate the different substructures (i.e., soft and semi-hard, or 2-jet and 3-jet events) of 
the distribution in the two classes of collisions, a will be the weight. 
The joint distribution for n F and n B charged particles is 

P**Ap>F, u b) = «Pi(n F , n B ) + (1 - a)P 2 (n F , n B ). (9) 

In the following, the term 'total' will be used used for quantities referring to the superposition 
of the two components, and the subscripts F and B stay as usual for 'Forward' and 'Backward' . 

Since we have defined the F and B hemispheres in a symmetric way, and since the collisions 
we are studying do not imply a difference between F and B hemispheres, it is clear that for each 
class i = 1, 2 of events, the joint MD's are symmetric in their arguments: 

Pi(n F ,n B ) = Pi(n B ,n F ). (10) 

in particular this implies that the average F multiplicity, np,i, equals the average B multiplicity, 
n Bji , and both are equal to half the average multiplicity fij in the two hemispheres, being rii = 

riF,i + n B>i : 

fiF,i = n B ,i = r^i- (H) 
We will also use below the equality of the variances: 

Dl Fji = (n F;i ) - n% = (n 2 Bji ) - n 2 B>i = (12) 
which together with Eq. (Q) applied to each class of events gives: 

{{n F ,i - n F ,i){n B>i - n B>i )) = hD 2 np i = biD 2 nB i ] (13) 
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the relation with the total variance then follows: 

Dli = Dl Fti + D 2 nBii + 2({n Fti - n F J{n B ,i - n B}i )) = 2(1 + k)D 2 nF ^ (14) 

Using Eq. (Q), we obtain the following relations for the total average multiplicity and vari- 
ance: 

n F = n B = -n; (15) 
Dl F = D 2 nB . (16) 

We start now by calculating the total forward average multiplicity n F : 

n F = an F> \ + (1 — a)n F;2 . (17) 

We proceed to calculate the total forward variance D 2 , using appropriately all the above men- 
tioned relations: 

Dl F = (n 2 F ) - nl 

= a(n 2 FA ) + (1 - a)(n^ 2 ) - [an F1 + (1 - a)n F2 ] 2 
= aD 2 np:1 + (1 - a)D 2 nF>2 + a{\ - a) (n F2 - n Fl f (18) 
aDl A (1 - a)D 2 2 1 2 

= 2[TTh) + 2(1 + **) • 

The total covariance can be calculated as follows: 

((n F - n F )(n B - n B )) = a(n F1 n BA ) + (1 - a){n F2 n B ^ 2 ) 

- [an F i + (1 - a)n F2 ] [an B}1 + (1 - a)n B>2 ] 

= a((n F1 - n F>1 )(n Bjl - n Bjl )} + (1 - a){(n Fj2 - n F>2 )(n Bt2 - n B>2 )) 

+ a(l - a) (n F2 - n F1 ) 2 (19) 

= abiD 2 np l + (1 - a)b 2 D 2 np2 + a(l - a) {n F2 - n F1 f 

aD 2 nl (l-a)D 2 n2 1 _ 2 

— F>7i rT~^2 + 7«(1 — a) (n 2 — ni) . 



2(1 + 6i) 2(1 + 63) 4 

The total correlation strength b for the weighted superposition of the two classes of events 
can therefore be written as: 

_ ahDl ^l + 6 2 ) + (1 - a)b 2 Dl 2 {l + &Q + fa(l - a){n 2 - n^jl + 6Q(1 + 6 2 ) 
~ 0^(1 + 62) + (1 - oyoyi + h) + ±a(l - a)(n 2 - rn) 2 (l + 6i)(l + b 2 ) 

(20) 

It should be remarked that even if the correlation strengths for events of class 1 and 2, bi and b 2 , 
separately vanish, the total strength does not: it depends on the difference in average multiplicity 
between the events of the two classes, namely 

b = \a{l - a){n 2 - fir) 2 

12 a/) 2 1 + (1 - ®) D l,2 + |a(! - oi)(fh -«i) 2 ' 
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Notice that in case each component is of NBMD type, the above formulae can be rewritten in 
terms of the standard NBMD parameters according to Eqs. (0) and @. 

We conclude that for b\ = 62 = forward-backward correlations are different from zero: 
the superposition of events of different classes generates a certain amount of positive FB corre- 
lations. 

It should be stressed that up to now all results have been obtained, as announced at the 
beginning of this Section, independently of any specific form of the charged particle MD's of 
the classes of events 1 and 2 contributing to the total charged particle MD and that Eq. ( [21] ) 
gives the amount of the superposition effect of the substructures 1 and 2, 6 12 , to the correlation 
strength with b\ = 62 = in terms of a, ni, D\ , n 2 , D\ parameters only. Once the different 
classes of events have been isolated and the above mentioned parameters measured, the estimate 
of 612 can easily be done. In e + e~ annihilation we are exactly in this situation and our test can 
be performed. It was shown by the OPAL Collaboration that for the 2-jet and the 3-jet event 
samples, separately, there is no correlations. Still there is correlation in the total sample, with 
b = 0.103 ± 0.007. Using a fit to OPAL data tf\X$ with similar conditions of the jet finder 
algorithm (a = 0.463, n x = 18.4, D 2 nl = 25.6, n 2 = 24.0, D 2 n2 = 44.6) we obtain from 
Eq. ( pit ) the value 6 12 = 0.101, in perfect agreement with the data. 

On the other hand, using UA5 two-components results [ [TT| ] (a = 0.75, fi\ = 24.0, D 2 1 = 
106, n 2 = 47.6, D 2 2 = 209) in full phase-space (thus including more correlations than present 
in the data, but nothing better is available in the actual phase-space range used by UA5) we 
obtain 6x2 = 0.28, a number much lower than the value found in the experiment (0.58). 

We conclude that FB correlations generated by the superposition of events of different 
classes are enough to explain observed FB correlations in e + e _ annihilation but not in hh 
collisions, where there exists a certain amount of correlation left within each class of events 
which should be taken into account. These results are a striking proof of the existence of the 
superposition effect, which was up to now only a guess, and of its relevance. 



3 Clan production and the correlation strength 

The next step in our approach is to calculate correlation strengths b\ and 6 2 in Eq. ( |2"0D in the two 
different classes of events in pp collisions in order to reproduce experimental data on b, which 
— as pointed out at the end of the previous Section — are not correctly reproduced by the knowl- 
edge of 612 only. The success of the superposition mechanism of weighted NB(Pascal)MD's for 
describing anomalies found in pp collisions strongly suggest to proceed to calculate b\ and b 2 
by using NB properties. 

In a naive approach to the problem one can try to apply Eq. (g), which gives at 546 GeV 
cm. energy b = 0.78: a much larger value than the experimental one (0.58). We conclude that 
charged particles FB distribution is not compatible with independent emission but is compatible 
with the production in clusters. Within the framework of the NBMD, these will be identified 
with clans. 

As already mentioned, clans can be produced forward or backward and may or may not leak 
particles to the opposite hemisphere. 
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3.1 The no-leakage case 

First, we will treat the general case in which no assumption is made about the MD of clans 
and of particles within a clan; then we will specialise our results to the case of the NBMD 
(Poissonian clans, logarithmic MD within a clan) obtaining a fair simplification of all formulae. 
This treatment is valid for each component separately, but the component index is dropped here 
to simplify the notation; Eq. ( |21]| ) can be used to obtain the total correlation strength. 

Accordingly, we write the joint distribution in n F ,n B as a convolution over the number of 
produced clans: 

P(n F ,n B )= V(N F ,N B )p F (n F \N F )p B (n B \N B ), (22) 

N F ,N B 

where we have indicate with capital N the number of clans and with V the joint distribution for 
N F clans forward and N B clans backward. Here p F (n\N) is the forward particle multiplicity 
distribution conditional on the number of forward clans, which, by arguments of symmetry, is 
the same distribution as p B (n\N). 

The symmetry of the reaction and of the hemispheres definition imply some conditions on 
V, namely that the average number of F and B clans at fixed full number of clans N are equal, 
and similarly for the corresponding variances: 

N F (N) = N B (N) = N/2; (23) 
d% F (N) = d% B (N) = (NUN)} - iV 2 /4. (24) 

If we now indicate with q(n) the MD within one clan, i.e., we write: 

PF(n\N) = ^2---^2q(n 1 )---q{n N y, (25) 

ni njv 
{Ei™i=™} 

then it is straightforward to show that the average value and the variance of n at fixed iV equal 
Nn c and ND 2 respectively. Using these results, it can be shown that 

(n F n B ) = (N F N B )n 2 c = Q(iV 2 ) - (d% F (N))^ n\ (26) 

and 

(n 2 F ) = (N F (N))Dl + (N F (N))n 2 c = \nD\ + n\ Ud\ p {N)) + \(N 2 )^ . (27) 

Then, without making any additional hypothesis on the clan distributions used, we can cal- 
culate the variances in the particle multiplicity: 



((n F - n F f) = ((n B - n B f) = ^ND 2 C + n\ 



\d 2 n + (d 2 NF {N)) 



(28) 



where N is the average full number of clans and D 2 N = (N 2 )—N 2 is the variance of the full clan 
multiplicity distribution; d% (N) is the variance in the distribution of the number of forward 
clans given that the total number of clans is N, and (d% (N)} is its average over iV; n c is the 
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full average number of particles per clan, and D 2 is the corresponding variance. The covariance 
is 



({n F - n F )(n B - n B )) = n\ 
The final result is thus: 



\d% - (d% F (N)) 



D 2 N - 4(d 2 Np (N)) 



D% + A(d% F (N))+2ND^ 
This result can be expressed in terms of quantities referring to the full MD: since 



(29) 



(30) 



D 2 n = D 2 c N + D 2 N n 2 c} (31) 



we obtain 

u D 2 Jn-D 2 /n c -An c (d% F (N))/N 



(32) 



Dl/n + D 2 /n c + 4n c (d% F (N)) /N' 

From the above formula, assuming binomially distributed clans, one can deduce Eq. (Q). 

The last formulae are of general validity, but now we specialise it to the NBMD case, where 
the overall clan multiplicity distribution is Poissonian and thus one finds D 2 N = N, and where 
clans do not talk to each other in this framework, thus the forward distribution at fixed total 
number of clans is binomial, (d 2 Np (N)) — N/4; putting these features together in Eq. ( |3~0| ) we 
obtain: 

N -N 

0. (33) 



N + N + 2ND 2 /n 2 c 

It follows the theorem: If particles are grouped in Poisson distributed clans; if each clan falls in 
the forward or backward hemisphere with the same probability independently of the other clans; 
if there is no leakage of particles from one hemisphere to the other: then no forward-backward 
correlation exist. 

We conclude that in order to have FB correlation in the framework of the clan interpretation 
of the NBMD it is necessary to allow clans to leak particles from one hemisphere to the other, 
i.e., clans must extend rather far in rapidity. 



3.2 The extension to the leakage case 

We return now to the general treatment. We start by writing the joint distribution as a con- 
volution over the number of produced clans and over the partitions of forward and backward 
produced particles among the clans: 

P(n F ,n B ) = V{N F ,N B ) ^ PF(m' F ,m' B \N F )p B (m%,m" B \N B ), (34) 

Np,Ns m' F +m' F =n F 
m' B +m'^=n B 

where Pi(m' F , m! B \Ni) is the joint probability of producing m F F-particles and m B B-particles 
from Ni clans in the i hemisphere (i = F, B). Notice that the summations are constrained to 
n F F-particles and n B B-particles, respectively. The usual symmetry argument imply that 

p F (n, m\N) = p B (m, n\N). (35) 
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It is now straightforward to show that it is still true that n F = Nn c /2. Similarly to Eq. (| 
we write now the decomposition 

p F (m F , m B \N F ) = 22 1f(^f,i, Hb,i) ■ ■ ■ <1f(hf,n f , Vb,n f ), (36) 

where q F ([i F , ^b) is the joint probability of producing, within one clan, \x F F-particles and hb 
B-particles. Of course one has: 

^2 QF(jiF,fiB)=q(n c ). (37) 



Define now 

pn c = ^2fi F q F (fj, F ,p, B ) = ^2,^ B q B {^ F ,^ B ), (38) 
qn c = ^HBqFiHF, I^b) = y^ j fJ, F q B {fJ, F , Hb), (39) 

where the second equality in each row has been obtained from the usual symmetry arguments, 
which also imply p + q = 1. The just defined parameter p controls the leakage from one 
hemisphere to the other: p = 1 means that no particle leaks, while 0.5 < p < 1 indicates 
leakage (p cannot be smaller than 0.5, or else the clan is classified in the wrong hemisphere). 

An important role in the final formula will be played by the covariance, 7, of fi F forward 
and fi B backward particles within a clan: 

7 = ((fi F - p F ) (/ab -Pb)) = ^2fj-Ffj-BqF(fj-F,fj-B) -pqn c - (40) 

In general, this quantity cannot be expressed in terms of n c or D 2 unless some explicit distribu- 
tion for the forward distribution at fixed number of particles per clan is assumed; for example, 
when particles within one clan are F-B binomially distributed, then 7 = [D\ — n c )pq. 
One then finds for the variance 



{{n F - npf) = ((n B - n B ) 2 ) = l -ND 2 c -N>y + n\ 



\D% + (d% F (N))(p-q) 2 



(41) 



-compare with Eq. ( j28[ ) — and for the covariance 

"1 



((n F - n F )(n B - n B )) = N>y + n 



^ N -{dj, F (N))(p- q y 



(42) 



-compare with Eq. 
The final general result for clans, for each component, is thus 

Dl-A^mtP-qf + ^/nl 



c 



D% + 4(d 2 NF (N))(p - q) 2 - ANj/n 2 c + 2ND 2 /n 
D 2 Jn - D 2 /n c - A(d 2 NF (N))(p - qfnjN + 4 7 /n e 
D 2 Jn + D 2 /n c + A(d 2 NF (N))(p - q) 2 n c /N - 4 7 /n c ' 



(43) 



9 



For NBMD clans, i.e., Poissonian (D 2 N = N) and independent (4(d 2 NF (N)) 
has inside a clan a logarithmic distribution, for which 



Dl 



-71, 



log(l 



N) clans, one 



(44) 



and 



V 



(V - 1) log(l - b') ' 



with 



(45) 



n 



n + k 



The correlation strength can then be written as 



{p-qf + 47/n; 



c 



(46) 



-l + (p-q) 2 - 4 7 /n2 - 2 log(l - b') jb' ' 



(47) 



As anticipated, the correlation strength is no longer zero. Notice that the no-leakage case can 
be re-obtained by setting p = 1 and consequently q = 0, 7 = 0. Recall that 7 is the covariance 
between the forward and backward multiplicities in one clan, so in Eq. ( p7| ) only parameters 
related to the within-clan distributions appear. If one were to assume that particles within a clan 
are independently distributed in the two hemispheres, so that 7 = [D 2 — n c )pq, then one arrives 
at a very simple formula: 



2b' pq — 1 + p 

l-2b'pq = 



(48) 



Notice that when p = 1/2 we recover the expression for the NBMD with binomially distributed 
particles, Eq. @. For b' — > 1, b turns out to depend on the p parameter only, a fact which will 
be very useful in the next Section. 



4 Energy dependence of the correlation strength 

In Fig. [1|, the behaviour of &i 2 as a function of cm. energy, as suggested by Eq. (plj), is shown 
in the framework of the three scenarios proposed by the present Authors in Ref. Jl2|]. The 
scenarios are based on the weighted superposition of soft (without mini-jets) and semi-hard 
(with mini-jets) events in pp collisions, each class of events being described by a NB(Pascal)MD 
with different characteristic parameters n and k. In the first scenario KNO scaling is assumed for 
both components; in scenario 2 KNO scaling is strongly violated for the semi-hard component 
and satisfied for the soft one; the third scenario is a QCD inspired scenario and its predictions 
turn out to be in general intermediate between the previous two. In addition, two alternative 
cm. energy dependences have been proposed for the semi-hard component: the first one is a 
consequence of the UA1 analysis on mini-jets and leads to 

ThWs) ~ 2ni(VS); (49A) 
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Figure 1: Predictions for the correlation coefficient for the superposition of two components in pp collisions in the 
case in which each component by itself presents no correlations; prediction are given in the three scenarios in full 
phase-space as function of the cm. energy. The dotted line is a fit to experimental values [Q]. 

the second one postulates that n 2 increases more rapidly with cm. energy (it takes into ac- 
count an eventual high particle density production in the central rapidity region) and correct the 
previous equation as follows 

n 2 {^) ~ 2ni(Vi) + c'ln 2 (v^); (49B) 



the estimate of d from existing fits is ps 0.1. The general trend of n as described by Eq. ( |49BD 
seems to be favoured in the present approach. 

It should be noticed that the superposition effect alone does not reproduce the logarithmic 
energy dependence of the correlation strength b in the GeV and TeV regions. 612 does not de- 
pend on pi and p 2 parameters, it is an increasing function of cm. energy in the GeV region 
hardly distinguishable in the three scenarios. At 900 GeV cm. energy, b 12 general trends over- 
lap, then they start to decrease smoothly and to differentiate their behaviour in the TeV region. 
Accordingly, in all scenarios Eq. (|20~1), and not Eq. (|T]), should be used in order to get the 
correct b behaviour. The role of b X2 has been shown to be fundamental in understanding FB 
correlations in e + e~ annihilation when particle population within each clan are quite small 
1-2) and particle leakage from one hemisphere to the other quite an exceptional fact (remember 
that here b\ ~ b 2 ~ 0). In proton-proton collisions the decrease of b\ 2 with cm. energy goes 
together with the onset of a much larger leakage activity from clans with large particle popu- 
lation, characteristic in particular of semi-hard events in scenarios 2 and 3. It is interesting to 
remark indeed in Fig. [T] that the decrease of b\ 2 in the TeV region is more pronounced in these 
two scenarios than in scenario 1 and in scenario 2 with respect to scenario 3, i.e., in general 
when clans with larger number of particles are produced and leakage effect is expected to be 
more important. 

Figure [2] summarises our findings on the cm. energy dependence of the FB charged par- 
ticle multiplicity correlation strength parameter b both for the soft and semi-hard components 
individually and for their superposition in the GeV and TeV regions up to 14 TeV for the three 
above mentioned scenarios. By assuming that in pp collisions in the ISR energy range the semi- 
hard component is negligible with respect to the soft one and by using the experimental b value 
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Figure 2: Predictions for the correlation coefficients for each component (soft and semi-hard) and for the total 
distribution in pp collisions. For each scenario, three cases are illustrated, corresponding to the three numbered 
branches: leakage increasing with yfs (upper branch, ©), constant leakage (middle branch, ©) and leakage de- 
creasing with ^/s (lower branch, ©). Leakage for the soft component is assumed constant at all energies. The 

dotted line is a fit to experimental values [H]. 



at 63 GeV in Eq. (]48|), the parameter p x controlling the leakage of particles emitted by clans in 
one hemisphere to the opposite one can be determined: we find p\ = 0.78, i.e., on the average 
22% of particle population within a clan are expected to leak in the opposite hemisphere. Since 
the average number of particles per clan for the soft component goes from « 2 at 63 GeV up to 
« 2.44 at 900 GeV (clans are almost of the same size) it is quite reasonable to conclude that p\ 
is approximately energy independent in the region. 

By using the p\ value for the soft component at 546 GeV in Eq. d48[), the value of the 
parameter p 2 controlling leakage effect in the semi-hard component can also be determined: 
it is found that p 2 = 0.77, a pretty close value to p\. Clan structure analysis can help in 
understanding this result: average clan size goes from 1.64 at 200 GeV up to 2.63 at 900 GeV 
for the semi-hard component. A very small difference with respect to the clan size for the soft 
component at the same cm. energy. In addition, since the increase of clan size for the semi- 
hard component from 200 GeV up to 900 GeV is indeed quite small, the p 2 value 0.77 can be 
considered approximately also constant in the GeV energy range. 

With the just mentioned constant values of p\ and p 2 parameters, and assuming the cm. 
energy dependence of the average charged particle multiplicity n corrected by a In 2 s term 



as indicated in Eq. Q49B[ ), the general trend of the b energy dependence from ISR up to top 



pp CERN Collider energy obtained by superimposing soft and semi-hard components effects is 
correctly reproduced. It agrees in particular with the phenomenological fit proposed by different 
Collaborations [§, |9|] in the full range, i.e., b = —0.019 + 0.061 In s (dotted line in Figs. [I] and 

The just mentioned assumptions on pi, p 2 and n are extended to the TeV region. At 900 GeV 
a clear bending in the behaviour of the total strength is visible in all three scenarios (full line © in 
Fig. 0) and is not compatible with the logarithmic increase of b in the GeV energy range. In view 
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of this result, should one question the validity of our assumptions on the constant behaviour of 
Pi and p 2 leakage parameters for the soft and semi-hard components in the TeV region? Lack of 
sound experimental data prevent us from making sharp statements on the problem. A possible 
insight comes again from clan structure analysis in the new energy domain. 

The average number of particle per clan for the soft component is growing from 2.63 at 900 
GeV up to 2.98 at 14 TeV (and from 2 at 63 GeV). Therefore constant p 1 seems a quite well 
founded assumption (it should be pointed out that variations of one-two per cent in the pi and 
P2 values as will be the case by increasing the average number of particles per clan of one unit 
do not change the overall scenario we are discussing). 

The same conclusion can be drawn in the TeV region for the semi-hard component in sce- 
nario 1 where the average number of particles per clans goes from 2.63 at 900 GeV to 3.28 at 
14 TeV and which shares with the soft component KNO scaling properties. But it is not true 
in general when KNO scaling violations are expected to occur: the average number of particle 
per clan goes for instance in scenario 2 from 2.63 at 900 GeV up to 7.36 at 14 TeV: one should 
see more particle leakage from clans and accordingly a decreasing p 2 value as cm. energy in- 
creases. This request can be taken into account as shown in Fig. ^: b bending becomes even 
higher assuming that p 2 is for instance a logarithmic increasing function of cm. energy in the 
TeV region and the corresponding leakage effect decreasing (curves ® in the figure). The op- 
posite situation occurs in the three scenarios and in particular in scenario 2 and 3 assuming that 
p 2 is a logarithmically decreasing function of the cm. energy (curves © in the figure) 

A similar global b behaviour is found indeed in the three scenarios when corresponding 
D^/ n 2 values are approximately the same. Larger differences appear when the increase of n 
with cm. energy is not compensated by a comparable decreasing of k parameter. It is also clear 
that asymptotically the above mentioned ratio will depend on k parameter only. 

The bending of the total FB multiplicity correlation strength in Fig. 0, in view of Eq. (^8j), 
is a natural consequence of the fact that the quantity n/(n + k) goes to 1 for increasing energy. 
This consideration notwithstanding, it is shown that an increasing leakage effect leads to an 
increase of b towards its maximum value 6=1. The impression is that in order to get a quick 
saturation of b, larger leakage and clans with higher population densities are needed: this could 
be a signal of a possible onset of a third class of hard events, harder than soft and semi-hard 
events discussed in this paper, producing the two requested effects. 

Although our approach is limited to full phase- space only and therefore kinematical con- 
straints at the border of the allowed rapidity range as well as short range correlations effects 
from central rapidity intervals might influence our claims, it seems quite reasonable to say that 
the bending effect should be expected in the TeV region. The deep connection between FB 
charged particle multiplicity correlation strength b, leakage effects and clan structure analysis 
in the framework of the superposition mechanism of different classes of events suggested by 
the present approach in order to try to understand and regulate b bending should also be tested 
in future experiment at LHC with the Alice detector. Our results on b bending effect should be 
compared with other predictions based on different models [ ]13| , [14]] • 

Another aspect of the problem should be examined for completeness: is the linear relations 
between Tlb^f) and rip, shown in Eq. (|I]), affected by the assumptions which have been made 
in order to determine the energy dependence of the total FB multiplicity correlation strength 6? 

Following the von Bahr-Ekspong theorem []T5|], in fact, the linear behaviour of ne(np) vs. 
riF, the binomial distribution of particles in the forward and backward hemispheres and the 
occurrence of the NBMD are not independent statements. They are strongly linked: the validity 



13 




Figure 3: Results of our model for nB{np) vs. np compared to experimental data ^[ in full phase-space at 63 
GeV (a) and in the pseudo-rapidity interval \rf\ < 4 at 900 GeV (b). Theoretical predictions at 14 TeV in f.p.s. for 
constant p2 (c) and p2 decreasing logarithmically with energy (d) are also shown. 



of any two of them implies the validity of the third. It is clear that the linearity of the relation 
is violated in the soft and semi-hard components separately in our approach: in fact, it should 
be pointed out that in the first step of the production process clans, not particles, are binomially 
distributed in the two hemispheres with p = 1/2, and that in the second step each clan generates 
particles, binomially distributed in the two hemispheres but with p x and p 2 different from 1/2. 
Accordingly, the final particles are not binomially distributed in the two hemispheres. Being the 
multiplicity distribution of the two classes of events a NBMD, the linearity of FB multiplicity 
correlations in each substructure cannot be exact. An example of a mild violation of linearity 
in the soft component is shown in Fig. ^|(a) where n^np) is plotted vs. n F at ISR energies 
(63 GeV). In addition, the extremely good fit of experimental data obtained by using a single 
NB(Pascal)MD for describing this component which in this case represents the total MD (at this 
energy the semi-hard component has been assumed to be negligible) is an indirect confirmation 
of the validity of our argument. Of course by going to higher cm. energies in pp collisions one 
expects that the linearity violation occur not only in the soft component but also in the semi-hard 
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one. The linearity violation in each separate components is shown in Fig. |3|(b) at 900 GeV cm 
energy: parameters from the two fits proposed in [ |TT| ] nave been used; their values determine 
the borders of the bands shown in the figure; the p 2 parameter has been taken equal to 0.77. 

It should be pointed out that since the total multiplicity distribution resulting from the 
weighted superposition of events of the two separate substructures (with and without mini-jets) 
is not of NB type, and in view of the lack of binomial structure in the final charged particles 
MD in the two hemispheres, the linear behaviour of fiB{n F ) vs. n F can in principle be restored 
for the total MD in agreement with the von Bahr-Ekspong theorem. 

In Fig. |](c) and (d), predictions at 14 TeV by using the parameters of the QCD-inspired 
scenario (labelled 3B in our notation) are given for p 2 = 0.77 as in the GeV region and for 
p 2 logarithmically decreasing with cm. energy in the TeV region (p 2 = 0.72 at 14 TeV). They 
confirm the general trend of n B {n F ) vs. n F seen at 900 GeV. The previously discussed scenarios 
IB and 2B lead in this case to very small modification (not shown). 

Some comments are needed in order to understand the theoretical predictions of the weighted 
superposition mechanism of the soft and semi-hard components for the general trend of (rip) 
vs.n F and its linear behaviour found experimentally in. (77] < 4 at 900 GeV by UA5 Collabora- 
tion (Fig. |(b)). 

The theoretical predictions are based on a separation of the two components which comes 
from fits [|TT]j to the total charged MD and not from an event-by-event classification. The idea 
of separating the total class of events into two components is correct in principle — as we have 
shown in our work — but is questionable here in its application. Different sets of NB parameters 
for the two components could lead to good fits of the total charged MD. Two sets we used 
are taken from [|TT]j and their predictions collected in bands in Fig. ^fb). Results should be 
considered therefore only indicative of the general trend, which we consider quite satisfactory, 
especially for large n F values in the semi-hard component and small n F values in the soft one, 
as visible in the figure. 

The second warning comes from intrinsic limitations of our approach. The choice to per- 
form our study analytically led us to consider the pi and p 2 parameters approximately constant 
throughout the GeV region, an approximation which might turn out to be not correct and in any 
case to imply small modifications of our results. 

In addition our predictions are based on the separation of the events into two classes. The 
onset of a third class of the same kind as that which could modify the b bending effect described 
in Fig. ^ — when properly inserted in our approach — could lead to different results. 

All these considerations apply of course also to the FB multiplicity correlation strength 
energy dependence discussed in Fig. || and testify that our predictions, in view of the lack of 
detailed experimental data, should be considerd only indicative of the expected general trends. 



5 Conclusions 

General formulas for FB multiplicity correlation strength in the framework of the superposition 
mechanism of two weighted MD's for different classes of events as functions of the average 
charged particle multiplicities and dispersions of each class of events are given. Assuming NB 
regularity behaviour for 2- and 3-jet samples of events in e + e~ annihilation at LEP energy, 
results obtained by OPAL collaborations are correctly reproduced within experimental errors in 
terms of the pure superposition effect, being FB multiplicity correlation strengths for the two 
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separate classes of events negligible. 

The same approach has been successfully extended to pp collisions. Differently from e + e~ 
annihilation, the FB multiplicity correlation strengths in the two substructures (soft and semi- 
hard events) turn out to be quite important and to lead to interesting predictions in the TeV region 
in the three scenarios discussed by the Authors in a previous paper and based on extrapolations 
of pp collisions properties in the GeV region. 

In particular an interesting connection is found between the particle populations within 
clans, particle leakage from clans in one hemisphere to the opposite hemisphere and super- 
position effect between different substructure of the collision. This finding favours structures 
with larger particle populations per clans and the decrease of the average number of clans. The 
FB charged particle multiplicity correlation strength is predicted to bend in all scenarios in the 
TeV region. 

The effects of the main assumptions of the present approach on the linear relation between 
the average number of particles emitted in one hemisphere as function of particle emitted in the 
opposite hemisphere have also been studied and their expected general trend at LHC explored. 
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